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Abstract 

We discuss ongoing progress towards precise characterizations of parton distributions at fixed 
transverse momentum, focusing on matrix elements in coordinate space and the treatment 
of endpoint singularities. 



Based on talks given at the HERA and the LHC Workshop (Hamburg, March 2007) 
and at the DIS07 Workshop (Munich, April 2007). 
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We discuss ongoing progress towards precise characterizations of parton distributions 
at fixed transverse momentum, focusing on matrix elements in coordinate space and 
the treatment of endpoint singularities. 

Parton distributions unintegrated in transverse momentum are naturally denned for small 
x via high-energy factorization pQ . This relates off-shell matrix elements with physical cross 
sections at x — ► 0, and gives a well-prescribed method to introduce unintegrated parton 
distributions in a gauge-invariant manner. 

The question of how to characterize gauge-invariantly a k± distribution over the whole 
phase space, on the other hand, is more difficult and not yet fully answered. Its relevance 
was already emphasized long ago in the context of Sudakov processes [2], jet physics [3], 
exclusive production [1], spin physics [5]. Although a complete framework is still missing, 
much work is currently underway on this subject, see e.g. [HI El El [HI EB El H2]- The 
discussion that follows focuses on aspects related to the gauge-invariant operator matrix 
elements and regularization methods for lightcone divergences. 

To ensure gauge invariance, the approach commonly used is to generalize the matrix 
elements that serve to define ordinary parton distributions to the case of field operators at 
non-lightcone distances [5j [13] . This leads one to consider the matrix element for the quark 
distribution (Fig. [1]) 

f(y) = (P|^)Fj(n) 7 +Vo(n)V>(0)|P) (1) 

with the quark fields ip evaluated at distance y = {0,y~,y±) for arbitrary y~ and y±, and 
the eikonal-line operators V given by 




Figure 1: Quark distribution function in the target of momentum p. 
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V y (n) = Vexp ( ig s / dr n^A^(y + t n) 



(2) 



where n is the direction of the eikonal line and A is the gauge field. 

However, while the use of Eq. |T]) does not pose major problems at tree level, it becomes 
more subtle at the level of radiative corrections. Part of the subtleties are associated with 
incomplete KLN cancellations that come from measuring kj_ in the initial state [21 El- These 
may appear as uncancelled divergences near the endpoints for certain lightconc momentum 
components [15] . Another set of issues are associated with the integration over all transverse 
momenta, and involve the relation of unintegrated parton distributions with the ordinary 
ones [TBI HZ1 HB] and the treatment of ultraviolet divergences. As observed in [19] for the 
case of the Sudakov form factor, the choice of a particular regularization method for the 
lightconc divergences also affects integrated distributions and ultraviolet subtractions. 

In |10] these effects are examined by an explicit calculation at one loop using techniques 
for the expansion of nonlocal operators. The answer for the coordinate-space matrix element 
is analyzed in powers of y 2 , separating logarithmic contributions from long distances and 
short distances, 



h(y) = asCp P + [ dv 

7T In 1 - v 



d Airp. 2 
2> ' p- 



(3) 



where ji is the dimensional-regularization scale and p is an infrared mass regulator. The 
lightconc singularity v — > 1 corresponds to the exclusive phase-space boundary x = 1. The 
singularity cancels for ordinary parton distributions (first term in the right hand side of 
Eq. l[5])) but it is present, even at d 4 and finite p, in subsequent terms, which contribute 
to the unintegrated parton distribution [TO] . This is then treated on the same footing as a 
physical correlation function, to be expanded in terms of the ordinary parton distributions 
with nontrivial, perturbatively calculable coefficient functions |16( I17j. 
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Figure 2: Cut-off regularization for the quark matrix element. 

Traditionally the effect of endpoint singularities is suppressed by the use of a cut-off. This 
is likely the case, for instance, in existing Monte-Carlo event generators that implement 
unintegrated parton distributions [20J [21] [22] [23] [24]. A cut-off is also implemented in 
treatments I2] ll21f25] based on regularizing the parton-distribution matrix element by taking 
the eikonal line n to be non-lightlike (Fig. [2]) , combined with evolution equations in the cut- 
off parameter r\ = (j> ■ n) 2 /n 2 [3] 126] - One-loop formulas in coordinate space corresponding 
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to the regularization method of Fig. [5] are given in [TO] , 
at fixed fc_L of order 



This method leads to a cut-off in x 

(4) 



However, cut-off regularization is not very well-suited for applications beyond the leading 
order. Furthermore, as the two lightcone limits y 2 — > and n? — > do not commute, a 
residual dependence on the regularization parameter r\ is left after integrating in k± the 
distribution defined with the cut-off. The relation with the standard operator product 
expansion is therefore not so transparent. 

An alternative approach is based on the subtractive regularization method [TU [27] . As 
explained in |f4j , in this approach the eikonal n is kept in lightlike direction but the singu- 
larities are canceled by multiplicative, gauge-invariant factors given by eikonal-line vacuum 
expectation values. The matrix element with subtraction factors is pictured in Fig. [3] where 
y = (0, y~,0±), and u is the direction of an auxiliary (non- lightlike) eikonal that provides a 
gauge-invariant regulator near x = 1 and cancels in the matrix element at y± = [lOj . The 
form of the counterterms is simple in coordinate space, where it can be given in terms of 
compact all-order expressions. 



The subtractive method is more systematic than the cut-off, and likely more suitable 
for using unintegrated parton distributions at subleading-log level. It can be useful for 
incorporating the unintegrated formulation in parton shower approaches [T7J [UJ HE] ■ Also, 
subleading accuracy is needed for matching large-x contributions with calculations at small 
x [HI El |H] and in the Sudakov region [301 EH]. 

The techniques discussed above will be instrumental to analyze factorization and evo- 
lution for parton distributions with increased precision [JJ [H] . The issue of soft gluon 
exchanges with spectator partons is revisited in [7] for hard pp collisions. A potential break- 
down of factorization at high order of perturbation theory is discussed (N 3 LO correction to 
dihadron production, with two soft and one collinear partons), which would be of interest 
to verify by calculation. Also, it will be interesting to investigate how the argument of [7] is 
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Figure 3: Matrix clement with subtractive regularization. 
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modified if one takes account of destructive-interference effects due to soft gluon coherence. 
A better understanding of these issues will help improve the present accuracy in current 
phcnomenological studies of the effects of partons' transverse momentum [24l [29l [30] . 
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